A model for non-singular black hole collapse and evaporation 
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We study the formation of a black hole and its subsequent evaporation in a model employing a 
minisuperspace approach to loop quantum gravity. In previous work the static solution was obtained 
and shown to be singularity-free. Here, we examine the more realistic dynamical case by generalizing 
the static case with help of the Vaidya metric. We track the formation and evolution of trapped 
surfaces during collapse and evaporation and examine the buildup of quantum gravitationally caused 
stress-energy preventing the formation of a singularity. 



o 

PL, 
(N 



i 



> 
m 

(N 
00 

(N 

On 
O 



X 



I. INTRODUCTION 

The formation of black hole singularities is an in- 
evitable consequence of General Relativity. As instances 
of infinite energy density and tidal forces, black holes 
have made headlines, inspired science fiction movies, and 
were studied in thousands of research articles. It adds 
to the fascination that we know today black holes are 
not just a mathematically possible solution to Einstein's 
field equations, but part of Nature. Since more than a 
decade now, we have good evidence that our Milky Way, 
as other galaxies, hosts many stellar black holes as well 
as a supermassive black hole in its center. 

From the perspective of quantum gravity, black holes 
are of interest because of the infinite curvature towards 
their center which signals a breakdown of General Rel- 
ativity. It is an area where effects of quantum gravity 
are strong, and it is generally expected that these effects 
prevent the formation of the singularity. Since the black 
hole emits particles in the process of Hawking radiation 
[H, the horizon radius decreases. In the standard case 
it approaches the singularity until both, the singularity 
and the horizon, vanish in the endpoint of evaporation 
0. However, if the singularity docs not exist, this sce- 
nario cannot be correct. Since the singularity plays a 
central role for the causal space-time diagram, its ab- 
sence in the presence of quantum gravitational effects 
has consequences for the entire global structure Q , and 
the removal of the singularity is essential for resolving the 
black hole information loss problem Q. To understand 
the dynamics of the gravitational and matter fields, it is 
then necessary to have a concrete model. 

It is thus promising that it has been shown in a sim- 
plified version of loop quantum gravity, known as loop 
quantum cosmology (LQC) Q, a resolution of singular- 
ities, the big bang as well as the black hole singularity 
, can be achieved. The regular static black hole met- 
ric was recently derived in Q, and studied more closely 
in [lfj. A resolution of the black hole singularity was 
also obtained in an effective, noncommutative approach 
to quantum grav ity 11 1 and in asymptotically safe quan- 
tum gravity [13J. In another work [12J, a 2-dimcnsional 
model was used to study the evaporation process in the 
absence of a singularity. 

Here, we will use a 4-dimensional model based on the 



static solution derived in [9| and generalize it to a dy- 
namical case which then allows us to examine the causal 
structure. This generalization holds to good accuracy in 
all realistic scenarios. This approach should be under- 
stood not as an exact solution to a problem that requires 
knowledge of a full theory of quantum gravity, but as a 
plausible model based on preliminary studies that allows 
us to investigate the general features of such regular black 
hole solutions. 

Non-singular black holes were considered alr eady by 
Bardeen in the late 60s and have a long histor y [l4| - l29| . 
We will here use a procedure similar to that in |24l|. The 
rest of the paper is organized as follows. We start in 
the next section by recalling the regular static metric 
we will be using. In section IIIII we generalize it to a 
collapse scenario and discuss its properties. In section 
IIVI we summarize the thermodynamical properties and, 
in section |Vj add the evaporation process and construct 
the complete causal diagram. The signature of the metric 
is (—,+,+,+) and we use the unit convention ft = c = 
Gn = 1. 



II. THE REGULAR 
SCHWARZSCHILD-METRIC 

Let us first summarize the regular black hole metric 
that we will be using. 

Loop Quantum Gravity (LQG) is a candidate theory of 
quantum gravity. It is obtained from the canonical quan- 
tization of the Einstein equations written in terms of the 
Ashtekar variables that is in terms of an su(2) 3- 
dimensional connection A and a triad E. The result [3lJ 
is that the basis states of LQG are closed graphs the edges 
of which are labelled by irreducible su(2) representations 
and the vertices by su(2) intertwiners. Physically, the 
edges represent quanta of area with area 7?p \J j(j + 1), 
where j is the representation label on the edge (a half- 
integer), lp is the Planck length, and 7 is a parameter 
of order 1 called the Immerzi parameter. Vertices of the 
graph represent quanta of 3- volume. The important ob- 
servation to make here is that area is quantized and the 
smallest quanta of area possible has area 

The regular black hole metric that we will be using is 
derived from a simplified model of LQG 0- To obtain 
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this simplified model wc make the following assumptions. 
First of all, the number of variables is reduced by assum- 
ing spherical symmetry. Then, instead of all possible 
closed graphs, a regular lattice with edge lengths <5i and 
82 is used. The solution is then obtained dynamically 
inside the homogeneous region (inside the horizon where 
space is homogeneous but not static). Analytically con- 
tinuing the solution outside the horizon one finds that 
one can reduce the two free parameters by imposing that 
the minimum area present in the solution corresponds 
to the minimum area of LQG. The one remaining un- 
known constant 5 is a parameter of the model determin- 
ing the strength of deviations from the classical theory, 
and would have to be constrained by experiment. With 
the plausible expectation that the quantum graviational 
corrections become relevant only when the curvature is 
in the Planckian regime, corresponding to 6 < 1, outside 
the horizon the solution is the Schwarzschild solution up 
to negligible Planck-scale corrections which allows us to 
believe the legitimacy of the analytical extension outside 
the horizon. 

This quantum gravitationally corrected Schwarzschild 
metric can be expressed in the form 



ds 2 



-G(r)dt 2 



dr 2 

W) 



H(r)dn 



(1) 



with dfl = dQ 2 + sin 2 



and 



G(r) = 



F(r) 



H(r) 



(r — r + )(r — r_)(r + r*y 



(r — r + )(r — r_)r 4 
(r + r*) 2 (r 4 +a 2 ) ' 

„2 



r 



(2) 



Here, r+ = 2m and r_ = 2mP 2 are the two horizons, and 
r * = \Z r + r - = 2mP. P is the polymeric function P = 
+ e 2 - + e 2 + 1), with e < 1 the product of 

the Immirzi parameter (7) and the polymeric parameter 
(6). With this, it is also P <C 1, such that r_ and 7% are 
very close to r = 0. The area is equal to A m ; n /87r, 
^4-min being the minimum area gap of LQG. 

Note that in the above metric, r is only asymptoti- 
cally the usual radial coordinate since gee is not just r 2 . 
This choice of coordinates however has the advantage of 
easily revealing the properties of this metric as we will 
see. But first, most importantly, in the limit r — > 00 the 
deviations from the Schwarzschild-solution are of order 
Me 2 /r, where M is the usual ADM-mass: 

9M 
r 



F(r) -> 1 - 
H(r) -> 



2AI 



(3) 



The ADM mass is the mass inferred by an observer at flat 
asymptotic infinity; it is determined solely by the metric 



at asymptotic infinity. The parameter m in the solution 
is related to the mass M by M = m(l + P) 2 . 

If one now makes the coordinate transformation R = 
ao/r with the rescaling t = ir 2 /ao, and simultaneously 
substitutes R± = ao/r T , i?* = ao/r« one finds that the 
metric in the new coordinates has the same form as in the 
old coordinates and thus exhibits a very compelling type 
of self-duality with dual radius r = ^/oq. Looking at the 
angular part of the metric, one sees that this dual radius 
corresponds to a minimal possible surface element. It is 
then also clear that in the limit r — > 0, corresponding 
to R — > 00, the solution does not have a singularity, but 
instead has another asymptotically flat Schwarzschild re- 
gion. 

The causal diagram for this metric, shown in Fig [IJ 
then has two horizons and two pairs of asymptotically flat 
regions, A, A' and B,B', as opposed to one such pair in 
the standard case. In the region enclosed by the horizons, 
space- and timelikeness is interchanged. The horizon at 
r_|_ is a future horizon for observers in the asymptotically 
flat B,B' region and a past horizon for observers inside 
the two horizons. Similarly, the r_ horizon is a future 
horizon for observers inside the two horizons but a past 
horizon for observes in A, A'. If one computes the time 
it takes for a particle to reach r = 0, one finds that it 
takes infinitely long [Icj • The diagram shown in Fig Q] is 
not analytically complete, but should be read as being 
continued on the dotted horizons at the bottom and top. 

The metric in Eq. ^ is a solution of a quantum gravi- 
tationally corrected set of equations which, in the absence 
of quantum corrections e, a — > 0, reproduce Einstein's 
field equations. However, due to these quantum correc- 
tions, the above metric is no longer a vacuum-solution 
to Einstein's field equations. Instead, if one computes 
the Einstein-tensor and sets it equal to a source term 
G>„ = 8ttT^, one obtains an effective quantum gravi- 
tational stress-energy-tensor T^ v . The exact expressions 
for the components of T are somewhat unsightly and can 
be found in the appendix. For our purposes it is here suf- 
ficient to note that the entries are not positive definite 
and violate the positive energy condition which is one of 
the assumptions for the singularity theorems. 



III. COLLAPSE 

We will proceed by combining the static metric with a 
radially ingoing null-dust, such that we obtain a dynam- 
ical space-time for a black hole formed from such dust. 
Usually described by the Vaidya metric [32| . we will in 
this scenario have corrections to the Vaidya metric that 
are negligible in the asymptotic region, but avoid the 
formation of a singularity in the strong-curvature region. 
The metric constructed this way in the following is not 
a strict solution of the minusuperspace LQC equations. 
However, as long as the null-dust does not already display 
strong quantum gravitational effects by its mass profile, 
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FIG. 2. G T V as a function of r for radially ingoing radiation 
and m'(v) = 1. The solid line depicts the classical case for 
e, ao — > 0. The long dashed line is for m(v) = 20, (r* > y^Oo) 
and the short dashed line is for m(v) = 5, (r» < t/oo). All 
quantities are in Planck units. 



FIG. 1. Penrose diagram of the regular static black hole solu- 
tion with two asymptotically flat regions. The both horizons, 
located at r+ and r_ , are marked in blue and red respectively. 

this solution should hold to good accuracjQ. 

We start by making a coordinate transformation and 
rewrite the static space-time in terms of the ingoing null- 
coord inate v. It is defined by the relation dv = dt + 
dr/^J F(r)G(r), which can be solved to obtain an explicit 
expression for v. The metric then takes the form 



ds 2 = -G(r)dv 2 + 2\ —-4 drdv + H(r)dfl . (4) 
V F ( r ) 

Now we allow the mass m in the static solution to depend 
on the advanced time, m — > m(v). Thereby, we will 
assume the mass is zero before an initial value v a and 
that the mass stops increasing at Vb- We can then, as 
before, use the Einstein equations G = 8?rT to obtain 
the effective quantum gravitational stress-energy tensor 



1 It has been claimed in |33| that, counterintuitively, quantum 
gravitational effects could become important already at the hori- 
zon when the collapse proceeds slowly. However, since we are 
considering null-dust, the collapse is as fast as can possibly be 
and these considerations do not apply. 



T. T v v and T r r do not change when m(v) is no longer 
constant. The transverse pressure T e s = however has 
an additional term 

^.a , , nn ^.a , \ Pr 2 m'(v) , . 

T e 9 (m(v)) = T e (m) - — — , 5 

2-K{r + 2m(v)F)* 

where ml = dm/dv. Because of the ingoing radiation, the 
stress-energy-tensor now also has an additional non-zero 
component, T r v , which describes radially ingoing energy 
flux 

2(1 + P) 2 r i (r*-a 2 ) )(r-r*(v))m>(v) 

(a 2 +r^(r + n(v)) 3 ' {) 

Notice that also in the dynamical case, trapping horizons 
still occur where g rr = F(r,v) vanishes [3J, |35|, so we 
can continue to use the notation from the static case 
just that r±(v) and r*(i>) are now functions of v. The 
r-dependence of this component is depicted in Fig [2] 

This metric reduces to the Vaidya solutions at large ra- 
dius, or for e — > 0, ao — > 0. However, in the usual Vaidya 
solutions, the ingoing radiation creates a central singular- 
ity. But as we see here, with the quantum gravitational 
correction, the center remains regular. 

We note that the ingoing energy flux has two zeros, one 
at r = r*(v) and one at r — y/ao, and is negative between 
these. What happens is that the quantum gravitational 
correction works against the ingoing flux by making a 
negative contribution until the effective flux has dropped 
to zero at whatever is larger, the horizon's geometric 
mean r* or the location of the dual radius r = -Joq. 
The flux then remains dominated by the quantum grav- 
itational effects, avoiding a collapse, until it has passed 
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r* and the dual radius where it quickly approaches what 
looks like an outgoing energy flux to the observer in the 
second asymptotic region. 

Since in the second asymptotic region A, A' the mass 
assigned to the white hole is inversely proportional to the 
ADM mass at r = oo, the white hole's mass must be de- 
creasing, consistent with the outgoing (or rather through- 
falling) energy flux. In this process, the past horizon will 
move towards smaller R or larger r, respectively. 



IV. THERMODYNAMICS 

Let us now briefly summarize the findings about the 
thcrmodynamical properties of this black hole solution, 
discussed in more detail in (Toj . 

Particle creation can take place at the horizons r + and 
r_ where there is high blueshift when tracing back light- 
rays. However, if the vacuum at I~ in the black hole's 
asymptotic region B, B' is empty of particles as usual, 
then there will be no flux from particle creation at r_ 
to I + in the second asymptotic region A, A'. This is a 
consequence of causality and energy conservation, which 
we can see as follows. 

Consider there was a particle creation at r_ resulting 
in a flux of Hawking radiation towards R = oo. The 
background is the time-reversed black hole situation but 
the flux is not time-reversed. This would mean a decrease 
of the white hole's mass for the observer at R = oo. How- 
ever, since our metric is geodesically complete, the parti- 
cles emitted at the white hole's horizon r_ can be traced 
back all the way to 1~ in the black hole's asymptotic re- 
gion B, B' . We recall that the white hole's mass for the 
observer in the A, A' region is inversely proportional to 
the black hole's mass and see that this particle creation 
at r_ would contribute to an increase of the black hole's 
mass corresponding to the decrease of the white hole's 
mass. Since there is particle emission also at the other 
horizon r + , we would have to add both fluxes to obtain 
the net mass change. 

However, we do as usual have a choice for the initial 
vacuum state at T~ and we will assume as normally that 
the vacuum in the black hole's asymptotic past is empty. 
From the above explanation we see now that this can only 
be the case if there is no particle flux from r_ to the white 
hole's asymptotic region 1 + . To achieve this, we have to 
chose the vacuum at T~ in the white hole's asymptotic 
region A, A 1 such that it contains a constant flux into 
the white hole with the effect that there is no outgoing 
particle flux created at r_ . This is the time- reversed situ- 
ation of an evaporating black hole with an empty ingoing 
vacuum. This situation is mathematically consistent be- 
cause particle production in the curved background only 
tells us the relation between the ingoing and outgoing 
vacuum states, but not the vacuum states themselves. 
We thus chose the vacuum state at X~ in the white hole's 
asymptotic region A, A' such that at r_ there is no ad- 



ditional outgoing flux created 0- 

Thus, the evaporation proceeds through the Hawking 
emission at r+, and the black hole's Bekcnstcin-Hawking 
temperature, given in terms of the surface gravity k by 
T BH = k/2tt, yields 



(2to) 3 (1-P 2 ) 
TBH{m) 4n[(2 m y + al] 



(7) 



This temperature coincides with the Hawking temper- 
ature in the limit of large masses but goes to zero for 
to -> 0. 

The luminosity can be estimated by use of the Stefan- 
Boltzmann law L{m) = aAH(m)T^ H (m), where (for a 
single massless field with two degrees of freedom) a = 
7r 2 /60, and A H (m) = 47r[(2m) 2 + a^/ (2to) 2 ] is the surface 
area of the horizon. Inserting the temperature, we obtain 



L(m) 



16m lu a(l -P 2 ) 4 



n*(a 2 



16 m 4 ) 3 



The mass loss of the black hole is given by —L(m) 
dm(v) 



dv 



-L[m(v)\ 



(8) 



(9) 



and we can integrate its inverse to obtain the mass func- 
tion m(v). The result of this integration with initial con- 
dition m(v = 0) = toq is 



5a% + 432a^TO 4 + 34560a 2 1 TO 8 - 61440to 12 )tt 3 

720to 9 (1 - P 2 ) 4 a 
5a 6 + 432a t m 4 + 34560a 2 m* - 61440toJ 2 )tt 3 

720to ) (1 - P 2 ) 4 a 



(10) 



In the limit m — > this expression becomes v ss 
ao7r 3 /(144m 9 (l — P 2 ) 4 a), and one thus concludes that 
the black hole needs an infinite amount of time to com- 
pletely evaporate. 



V. COLLAPSE AND EVAPORATION 

We are now well prepared to combine formation and 
evaporation of the black hole. As in section llHl we divide 
space-time into regions of advanced time. We start with 
empty space before v a , let the mass increase from v a to Vb, 
and stop the increase thereafter. Hawking radiation will 
set in, but for astrophysical black holes this evaporation 
will proceed very slowly, such that we have a long time 
span during which the black hole is quasi-stable and m 
remains constant to good accuracy at too. Then, at some 
later time, v c , Hawking radiation becomes relevant and 



2 Alternatively, we could demand the vacuum at past infinity in 
the second asymptotic region to be free of particles, but then 
the vacuum in the black hole region's past infinity would have 
to contain particles. We will not further consider this possibility 
here. 
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to decreases until it reaches zero again. As we have seen 
in the previous section, it will reach zero only in the limit 
v —> oo. 

We thus have the partition — oo < v a < Vb < v c < oo 
with 

Vw G (—00, v a ) : m(v) = 0, (11) 
\fv G (v a , v b ) : d/dv m(v) > 0, (12) 
Vu G (v b ,v c ) : m(v) = m , (13) 

Vw G (u c , +00) : d/dv m(v) < 0, (14) 
for v — > +00 : m(v) — > 0. (15) 

Strictly speaking the mass would immediately start drop- 
ping without incoming energy flux and thus v a = Vb, but 
stretching this region out will be more illuminating to 
clearly depict the long time during which the hole is qua- 
sistable. 

To describe the Hawking-radiation we will consider the 
creation of (massless) particles on the horizon such that 
locally energy is conserved. We then have an ingoing 
radiation with negative energy balanced by outgoing ra- 
diation of positive energy. Both fluxes originate at the 
horizon and have the same mass profile which is given by 
the Hawking temperature. The area with ingoing neg- 
ative density is again described by an ingoing Vaidya 
solution, while the one with outgoing positive density is 
described by an outgoing Vaidya solution. 

The outgoing Vaidya solution has a mass-profile that 
depends on the retarded time u instead of v and the 
mass decreases instead of increases. The retarded time is 
defined by du — dt — drj \J 'F(r)G(r). After a coordinate 
transformation, the metric reads 




FIG. 3. Penrose diagram for the formation and evaporation 
of the regular black hole metric. The red and dark blue solid 
lines depict the two trapping horizons r_ and r+ . The brown, 
dotted line is the curve of r = y/ao and the black, short dashed 
one is r* . The light blue arrows represent positive energy flux, 
the magenta arrows negative energy flux. 



ds 2 



-G(r, u)du 2 



<G(r, u) 
F{r,u) 



dudr + H(r)dn , (16) 



where F(r, u) and G(r, u) have the same form as in the 
static case ©, but with where m is replaced by a func- 
tion m(u). We fix the zero point of the retarded time u 
so that r = r + corresponds to u c = v c . Then there is 
a static region with total mass mo for v > v c , u < u c . 
Note that since the spacetime described here has neither 
a singularity nor an event horizon, we can consider pair 
creation to happen directly at the trapping horizon in- 
stead of at a different timelike hypersurface outside the 
horizon, as done in [3|| . We have in this way further par- 
titioned spacetime in regions, broken down by retarded 
time: 



Vit < u c : m(u) = mo , 
Vu > u c '- d/du m(u) < 



(17) 
(18) 



Now that we have all parts together, let us explain the 
complete dynamics as depicted in the resulting causal 
diagram FigO 

In the region v < v a we have a flat and empty region, 
described by a piece of Minkowski-space. For all times 
v > v a , the inner and outer trapping horizons are present. 



These horizons join smoothly at r = in an infinite time 
and enclose a non-compact region of trapped surfaces. 

A black hole begins to form at v = v a from null dust 
which has collapsed completely at v = Vb to a static state 
with mass Too- It begins to evaporate at v = v c , and the 
complete evaporation takes an infinite amount of time. 
The observer at I + sees particle emission set in at some 
retarded time u c which corresponds to the lightlike sur- 
face where the horizon has lingered for a long time. The 
region with v > v c is then divided into a static region 
for u < u c , and the dynamic Vaidya region for u > u c , 
which is further subdivided into an ingoing and an out- 
going part. 

As previously mentioned, the radially ingoing flux 
(light blue arrows) in the collapse region is not positive 
everywhere due to the quantum gravitational contribu- 
tion. It has a flipped sign in the area between r* (black 
short dashed curve) and r = ^Joq (brown dotted curve) 
which is grey shaded in the figure. Likewise, the ingoing 
negative flux during evaporation (magenta arrows) has 
another such region with flipped sign. It is in this region, 
between the two horizon's geometric mean value and 
the dual radius corresponding to the minimal area, that 
the quantum gravitational corrections noticeably modify 
the classical and semi-classical case, first by preventing 
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the formation of a singularity, and then by decreasing the 
black hole's temperature towards zero. 

VI. CONCLUSIONS 

We have investigated a model for collapse and evapora- 
tion of a black hole that is entirely singularity-free. The 
spacetime docs not have an event horizon, but two trap- 
ping horizons. By generalizing the previously derived 
static metric to a dynamical one by use of the Vaidya 
metric we found that the gravitational stress-energy ten- 
sor builds up a negative contribution that violates the 
positive energy condition and prevents the formation of 
a singularity. We divided spacetime into six different re- 
gions described by different metrics, and constructed the 
causal diagram for the complete evaporation. The value 
of the scenario studied here is that it provides a concrete, 
calculable, model for how quantum gravitational effects 
alter the black hole spacetime. 
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coordinate (v,r,6,ip) are: 



v (a 2 + r A ) 3 (r + n) 3 
[2agr 4 (r + r*)(6r 2 + 7r_r + — 7r(r_ + r + ) — 2rr* 
~rl) ~ a^-r-r+r* + + 2r 2 (r_ + r+ + 2r„) 
+3r(-r_r+ + r 2 )) - r 6 {-r^r + i\ + rl 
+r(r-r+ + 2(r_ + r+)r* + 3r 2 ))] , 



(a2+r 4 )3(r + rv)3 
[2o^r 4 (r + r*)(-r_r + + r\ + r(r_ + r+ + 2r*)) 
+a 4 l (4r 3 - 2r 2 (r_ + r + - 2r*) - r_r+r* + r 3 
+r(r-r + + 3r 2 )) + r 8 (4r 2 r* + 3r_r+r* + r 3 
+r(r_r + - 2(r_ + r + )r» + 3r 2 ))] , 



6 (a 2 +r 4 )3( r + n )4 

[r 7 (r 2 r_r + + r(2r 2 + 6r_r + — 3r(r_ + r + ))r„ 

+ (r - 2r_)(r - 2r+)r 2 ) + 2a 2 r 4 (r 2 (r_ + r+) 

+ (r_ + r + )r 2 + r(r_ - r*)(— r + + r*)) 
+ao(4r 3 - 2r 2 (r_ + r + - 3r*) + 2r_r + r st 

+r(r_r+ - 3(r_ + r + )r* + rl))] . (19) 



APPENDIX 

The effective energy momentum tensor is defined by 
= G^/Stt. The components of the Einstein tensor in 
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